A directional specific heat study on the gap structure of overdoped Ba(Fei_ :c Co :r )2As2 
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Low-temperature specific heat is measured on the overdoped Ba(Fei_ a; Co :E )2As2 (x = 0.13) single 
crystal under magnetic fields along different directions. A clear anisotropy is observed on the 
electronic specific heat coefficient 'y(H). The value of j(H) is clearly larger with magnetic field 
along [001] (c-axis) than that within the a6-plane of the crystal lattice, which cannot be attributed 
to the effect by anisotropy of the upper critical field. Meanwhile, the data show a rather small 
difference when the direction of the field is rotated from [100] to [110] direction within the afc-plane. 
Our results suggest that a considerable part of the line nodes is not excited to contribute to the 
quasiparticle density of states by the field when the field is within the afe-plane. The constraint on 
the topology of the gap nodes is discussed based on the data. 

PACS numbers: 74.20. Rp, 74.70.Xa, 74.62.Dh, 65.40.Ba 



The gap structure of the Fe-based superconductors 
seems to be more complicated than expected, which 
is found to vary substantially from family to family. 
The consensus has been reached on several systems, 
e.g. LaFePO,^ KFe 2 As 2 ,^ BaFe 2 (Asi_ :c P a: )2r' S and 
so on, that nodes exist on the gap structure. How- 
ever, experimental results gave rather contradicting con- 
clusions on this issue in other systems of the Fe-based 
superconductors^— Furthermore, recent experiments 
reveal that the gap structure can modulate with the 
doping concentration even in the same family. In the 
case of electron-doped (Co- or Ni- doped) 122 system, 
the increase in the gap anisotropy and even the emer- 
gence of gap nodes, as the doping content increases in 
the overdoped region, are evidenced by many experimen- 
tal methods J^r— This behavior is attributed to the en- 
hancement of intraband interaction and pair scattering 
between the electron-like Fermi surfaces as the system is 
doped away from the optimal point. 22 

Specific heat (SH) is one of the most powerful bulk 
probes for investigating the gap structure of the un- 
conventional superconductors. The variation of elec- 
tronic SH (C' e i) versus temperature (T) and magnetic 
field (H) strongly depends on the gap structures of 
the superconducting materials. T dependence of C e i 
for a superconductor with line nodes varies as C e ; oc 
T 2 in low temperature ; 23 ' 24 while an exponential T 
dependence is expected for a superconductor with an 
isotropic gap. In the mixed state, H can induce vor- 
tices with a supercurrent flowing perpendicularly to the 
applied H. The low-energy quasiparticles will undergo 
a Doppler shift induced by the supercurrent, which can 
give rise to a considerable enhancement of C' e i for a nodal 
supcrconductor i 25 ' 26 This is called the Volovik effect. We 
note here that the Volovik effect is the strongest (weak- 
est) when H is parallel (perpendicular) to the Fermi sur- 



face where the nodes reside, because the Doppler shift 
energy is determined by 8E oc v s ■ Vf, where v s and vp 
are the supercurrent velocity and Fermi velocity respec- 
tively, and the Volovik effect works mainly at the places 
near the nodes. This provides valuable information on 
the topology of the line nodes, when we compare the SH 
data obtained with H along different directions of the 
crystal lattice. 

A clear evidence for the line nodes in the energy gap 
has been observed in terms of the T 2 term in C e i for the 
overdoped Ba(Fei_ :E Co :E )2As2 in our previous studies^ 
A quick increase of the electronic SH coefficient A-f(H) 
with H was observed when H was applied in the aft- 
plane, being consistent with the prediction from the 
Volovik effect on a nodal superconductor. Meanwhile, 
the underdoped and optimal doped samples show a small 
anisotropy on the energy gap. In order to obtain fur- 
ther information on the topology of the line nodes in the 
overdoped samples, directional SH studies are very war- 
ranted. Here we report such important experiments. The 
unambiguous anisotropic behaviors related to the Volovik 
effect imposed on line nodes can be confirmed by the fact 
that the electronic SH coefficient 7 (if) increases more 
quickly when if is parallel to the c-axis than that within 
the aft-plane. The possible constraint on the topology of 
line nodes is discussed based on the observations. 

The Ba(Fei_ a; Co a ;)2As2 single crystals were grown by a 
self-flux method. As-grown samples were annealed under 
high vacuum at 800 °C for 20 days. The dc magnetiza- 
tion measurements were made with a superconducting 
quantum interference device (Quantum Design, MPMS). 
SH was measured with if along three directions (see in- 
set (b) of Fig. 1). The data obtained with if within 
aft-plane were collected with a Helium-3 system attached 
with the physical property measurement system (Quan- 
tum Design, PPMS), while those with if along c-axis 
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FIG. 1: (color online) Main frame: temperature dependence 
of SH coefficient (C/T) for the overdoped Ba(Fei_ :E Co a; )2As2 
with x — 0.13 under zero field. Inset (a): diamagnetic transi- 
tion of the same sample. Inset (b): schematic representation 
of three directions along which the magnetic field is applied. 



were obtained in a standard mode. 

The superconducting transition of the selected sample 
with x = 0.13 is checked by the dc magnetization and 
specific heat measurements. The present sample was 
determined to be in the overdoped region of the phase 
diagram. 20 In Fig. 1, we show the T dependence of the 
SH coefficient C/T up to 15 K under zero field. The kink 
at about 12.5 K indicated by red arrow is the supercon- 
ducting transition. This temperature corresponds to the 
end point of the diamagnetic transition, as shown in the 
inset (a) of Fig. 1. 

We measured SH of Ba(Fc 1 _ ;r Co x )2As2 with x = 0.13 
under H along the [100], [110], and [001] directions of 
the crystal lattice. A schematic representation of the 
three directions is given in inset (b) of Fig. 1. The raw 
data of SH under H of these three directions are plotted 
as C/T vs T 2 in Figs. 2(a)-(c), respectively. Here we 
focus on the behaviors of our data in the low-T region 
below 4.5 K. One can see that the behaviors of the data 
with H along [100] and [110] directions are quite similar 
with each other. While the data with H along the [001] 
direction increase more quickly with H. Nevertheless, 
all the three sets of data show clear negative curvatures 
in the T range we studied. In our previous work^ we 
have attributed this behavior to the presence of the T 2 
term in the electronic SH, which is consistent with the 
prediction for the superconductors with line nodes in the 
energy gap. The data were then fitted by the following 
equation 



C(T, H) = j(H)T + a(H)T 2 + BT 3 



(1) 



where j(H) is the electronic SH coefficient under H, 
a(H) is the coefficient of the T 2 term under H, and 8 is 
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FIG. 2: (color online) Raw SH data of the Ba(Fei_ a: Co a ;)2As2 
sample with x — 0.13 under fields aligned with the [100], [110], 
and [001] directions of the crystal lattice, respectively. The 
dashed lines represent the results of theoretical fitting. 



the phonon SH coefficient. The value of B was found to be 
almost independent of H . Consequently, here we average 
the value of B under different fields and fix it when fitting 
the data using Eq. (1). The effect of small fluctuations 
of B is transferred to the error bar of the resulting fitting 
parameters ~f(H) and a(H) (see Fig. 3). The fitting re- 
sults are displayed by the dashed lines in Fig. 2. Only 
four selected fitting curves are shown in Figs. 2(a) and 
(b) respectively for clarity. It is clear that these curves 
describe the negative-curvature features commendably. 

The fitting parameters "f{H) and ct(H) are shown in 
Figs. 3(a) and (b). The horizontal ordinate is normalized 
with the upper critical field H C 2, so as to eliminate the 
effect of anisotropy of H C 2 on the H dependent data. As 
we have stated^ fl" c2 within the a6-plane is about 28 
T, from which the H C 2 value along the c-axis can be esti- 
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FIG. 3: (color online) Field dependence of y(H) and a(H) 
with H along three directions. H is normalized with the up- 
per critical field H C 2- The blue dashed line in (a) shows the 
theoretical curve for the case with an isotropic gap. 



mated using the its anisotropy (T = H^/H^ 2 ). Fr om the 
reported high-i/ experiments, we know that T decreases 
when reducing T and finally reaches about 1.1 at 0.7 K 
for the near-optimal doped Ba(Fei_ a Co a ;)2As2i 2 ''' 28 Here 
we take the value T ~ 1.3 for the present overdoped case, 
which would be no less than the actual value. Such a 
treatment will not affect the conclusions described below, 
because a smaller T would result in the enhancement of 
the anisotropic features. 

Both 7 (If) and a(H) show small deviations among the 
three directions when H is approaching zero, which con- 
firms the reliability of our data even though the data with 
H along c-axis is only measured down to 1.9 K. From 
Fig. 3(a) one can see that all the three sets of data are 
clearly above the blue dashed line, which represents the 
theoretical curve for the case with an isotropic gap. This 
is consistent with the features induced by the Volovik ef- 
fect in a nodal superconductor. It is clear that j(H) with 
H aligned with c-axis (hereafter abbreviated as 7 c (if)) 
increases more quickly than that within the afr-plane (ab- 
breviated as j ab (H)). We found that j ab (H) takes up 
about 70% of 7 c (-ff) when the reduced H/H C 2 equals to 
0.3. In sharp contrast, the data within the a6-plane re- 
main almost unchanged within the extent of error bar 
when H is rotated from [100] to [110] direction. Con- 
sidering the fact that the increase of j(H) is induced by 
Volovik effect, our present observations suggest that the 
Volovik effect is stronger when H is applied along the 
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FIG. 4: (color online) Schematic diagram of two typical mod- 
els of line nodes on one of the Fermi surfaces with a 3D charac- 
ter. Here we only display a longitudinal section of the Fermi 
surface shown by the two black lines. The vertical (a) and 
horizontal (b) line nodes are shown by the red lines. 



c-axis than that within the afe-plane. The anisotropic 
features of H dependent a(H) shown in Fig. 3(b) also 
support such an argument because the decrease of a(H) 
is associated with the Volovik effect. Obviously, these 
findings will provide a hint for investigating the topology 
of the nodes, which will be discuss in the following. 

Different topologies of the nodes or gap minima 
have been proposed and the significant influence by 
the three-dimensional (3D) features of the Fermi sur- 
face (FS) has been noticed, by both experiments and 
calculations ! 17 ' 18 ' 29 ' 30 Consequently, here we consider the 
constraint by our data on the line nodes based on the two 
typical models, where the line nodes with two different 
orientations are proposed on one of the FSs with a 3D 
character. As shown in Fig. 4, the longitudinal section 
of the FS is shown by two black lines. The red line in 
Fig. 4(a) and red circle in Fig. 4(b) represent the posi- 
tions of vertical and horizontal line nodes, respectively. 

We first check the situation of the vertical-node case 
as shown in Fig. 4(a). When H is applied parallel to 
k Zl which is shown by the blue arrow H-y, the segments 
of line nodes on the FS with no or small 3D features 
will experience a strong Volovik effect, while that with 
clear 3D features only undergo a depressed Volovik effect 
induced by the projection of H on the tangent surface 
(i?H = Hi cos#), where angle 9 describes the deviation of 
the FS from the direction of k z . When H is perpendicular 
to fc z , the situation will be rather complex because the 
direction can be rotated by 360°. The Volovik effect is 
the weakest if H is rotated to the nodal direction (see i?2 
in Fig. 4(a)) because only the segments of nodes on the 
3D-dispersed FS can experience the Volovik effect with 
a projected field H\\ — H2Sm9. The situations with H 
along other directions may vary depending on the number 
of the vertical line nodes and their distribution on the 
FS. Nevertheless, the clear anisotropy between j c (H) and 
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j ab (H) requires that the angle 9 or the proportion of the 
3D-dispersed FS should not be too large. 

One problem in the scenario described above is that 
"f ab {H) should show some variation in principle when H 
is rotated within a6-plane, which was not observed in our 
data. One possibility is that the number of nodes is large 
(e.g. lager than 4) , so that the variation will be weakened. 
Another explanation may be given by calculations based 
on an extended-s-wave casej^I where an attenuation of 
the SH vibration is predicted because of the elliptical FS 
pockets near the M points. Although both cases do not 
support the d-wave symmetry of the energy gap, we can't 
rule out the possibility that both [100] and [110] direc- 
tions deviate from the nodal direction, and consequently 
we failed to observe the difference between the two di- 
rections. This may need further clarification by detailed 
angle-resolved SH measurements. 

The situation becomes somewhat different if we have 
the horizontal line nodes, as shown in Fig. 4(b). When H 
is parallel to k z (see H3), the whole nodal line will experi- 
ence a depressed Volovik effect induced by H11 = i? 3 cos 9. 
Whereas when H is applied perpendicular to k z (see H4), 
the segments of nodes on the areas of the FS marked 
by the green circles will experience a strong Volovik ef- 
fect because H is roughly parallel to these parts of FS. 
Meanwhile, the segments of nodes marked by the violet 
circles will experience a depressed Volovik effect induced 
by £/]| = i/4sin6>. The fact that 7 c (£f) is clearly larger 
than 7 ab (if) in our data means that the overall Volovik 
effect when H is parallel to k z should exceed that with 
H perpendicular to k z . This implies that the angle 9 
cannot be too large. The advantage of this latter model 
is that the unobservable vibration of r y ab (H) can be ex- 



plained naturally. A recent angle resolved photoemission 
spectroscopy (ARPES) study on BaFe 2 (As .7Po.3)2 has 
revealed such a horizontal circular line node on the hole 
FS around the Z point at the Brillouin zone boundary. 32 
In summary, we studied the low-temperature SH on the 
overdoped Ba(Fei_ a; Co a ;)2As2 with H along three differ- 
ent directions. Clear anisotropic behaviors are observed 
from H dependent data. The electronic SH coefficient 
j(H) increases more quickly when H is in the c-axis di- 
rection than that within the a&-plane, whereas the data 
remain unchanged within our resolution when H is ro- 
tated within the afr-plane. Our results suggest that a 
considerable portion of the line nodes is not excited ple- 
narily to contribute to the density of states when H is in 
the ab in-plane. These conclusions supply an important 
constraint when investigating the topologies of the line 
nodes in this system. 
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